7 research outputs found
RED-PSM: Regularization by Denoising of Partially Separable Models for Dynamic Imaging
Dynamic imaging addresses the recovery of a time-varying 2D or 3D object at
each time instant using its undersampled measurements. In particular, in the
case of dynamic tomography, only a single projection at a single view angle may
be available at a time, making the problem severely ill-posed. In this work, we
propose an approach, RED-PSM, which combines for the first time two powerful
techniques to address this challenging imaging problem. The first, are
partially separable models, which have been used to efficiently introduce a
low-rank prior for the spatio-temporal object. The second is the recent
Regularization by Denoising (RED), which provides a flexible framework to
exploit the impressive performance of state-of-the-art image denoising
algorithms, for various inverse problems. We propose a partially separable
objective with RED and a computationally efficient and scalable optimization
scheme with variable splitting and ADMM. Theoretical analysis proves the
convergence of our objective to a value corresponding to a stationary point
satisfying the first-order optimality conditions. Convergence is accelerated by
a particular projection-domain-based initialization. We demonstrate the
performance and computational improvements of our proposed RED-PSM with a
learned image denoiser by comparing it to a recent deep-prior-based method
known as TD-DIP. Although the main focus is on dynamic tomography, we also show
the performance advantages of RED-PSM in a cardiac dynamic MRI setting
Machine Learning technique for isotopic determination of radioisotopes using HPGe -ray spectra
-ray spectroscopy is a quantitative, non-destructive
technique that may be utilized for the identification and quantitative isotopic
estimation of radionuclides. Traditional methods of isotopic determination have
various challenges that contribute to statistical and systematic uncertainties
in the estimated isotopics. Furthermore, these methods typically require
numerous pre-processing steps, and have only been rigorously tested in
laboratory settings with limited shielding. In this work, we examine the
application of a number of machine learning based regression algorithms as
alternatives to conventional approaches for analyzing -ray
spectroscopy data in the Emergency Response arena. This approach not only
eliminates many steps in the analysis procedure, and therefore offers potential
to reduce this source of systematic uncertainty, but is also shown to offer
comparable performance to conventional approaches in the Emergency Response
Application
Solving 3D Inverse Problems using Pre-trained 2D Diffusion Models
Diffusion models have emerged as the new state-of-the-art generative model
with high quality samples, with intriguing properties such as mode coverage and
high flexibility. They have also been shown to be effective inverse problem
solvers, acting as the prior of the distribution, while the information of the
forward model can be granted at the sampling stage. Nonetheless, as the
generative process remains in the same high dimensional (i.e. identical to data
dimension) space, the models have not been extended to 3D inverse problems due
to the extremely high memory and computational cost. In this paper, we combine
the ideas from the conventional model-based iterative reconstruction with the
modern diffusion models, which leads to a highly effective method for solving
3D medical image reconstruction tasks such as sparse-view tomography, limited
angle tomography, compressed sensing MRI from pre-trained 2D diffusion models.
In essence, we propose to augment the 2D diffusion prior with a model-based
prior in the remaining direction at test time, such that one can achieve
coherent reconstructions across all dimensions. Our method can be run in a
single commodity GPU, and establishes the new state-of-the-art, showing that
the proposed method can perform reconstructions of high fidelity and accuracy
even in the most extreme cases (e.g. 2-view 3D tomography). We further reveal
that the generalization capacity of the proposed method is surprisingly high,
and can be used to reconstruct volumes that are entirely different from the
training dataset.Comment: 14 pages, 10 figure
High-Precision Inversion of Dynamic Radiography Using Hydrodynamic Features
Radiography is often used to probe complex, evolving density fields in
dynamic systems and in so doing gain insight into the underlying physics. This
technique has been used in numerous fields including materials science, shock
physics, inertial confinement fusion, and other national security applications.
In many of these applications, however, complications resulting from noise,
scatter, complex beam dynamics, etc. prevent the reconstruction of density from
being accurate enough to identify the underlying physics with sufficient
confidence. As such, density reconstruction from static/dynamic radiography has
typically been limited to identifying discontinuous features such as cracks and
voids in a number of these applications.
In this work, we propose a fundamentally new approach to reconstructing
density from a temporal sequence of radiographic images. Using only the robust
features identifiable in radiographs, we combine them with the underlying
hydrodynamic equations of motion using a machine learning approach, namely,
conditional generative adversarial networks (cGAN), to determine the density
fields from a dynamic sequence of radiographs. Next, we seek to further enhance
the hydrodynamic consistency of the ML-based density reconstruction through a
process of parameter estimation and projection onto a hydrodynamic manifold. In
this context, we note that the distance from the hydrodynamic manifold given by
the training data to the test data in the parameter space considered both
serves as a diagnostic of the robustness of the predictions and serves to
augment the training database, with the expectation that the latter will
further reduce future density reconstruction errors. Finally, we demonstrate
the ability of this method to outperform a traditional radiographic
reconstruction in capturing allowable hydrodynamic paths even when relatively
small amounts of scatter are present.Comment: Submitted to Optics Expres
Ad hoc file systems for high-performance computing
Storage backends of parallel compute clusters are still based mostly on magnetic disks, while newer and faster storage technologies such as flash-based SSDs or non-volatile random access memory (NVRAM) are deployed within compute nodes. Including these new storage technologies into scientific workflows is unfortunately today a mostly manual task, and most scientists therefore do not take advantage of the faster storage media. One approach to systematically include nodelocal SSDs or NVRAMs into scientific workflows is to deploy ad hoc file systems over a set of compute nodes, which serve as temporary storage systems for single applications or longer-running campaigns. This paper presents results from the Dagstuhl Seminar 17202 “Challenges and Opportunities of User-Level File Systems for HPC” and discusses application scenarios as well as design strategies for ad hoc file systems using node-local storage media. The discussion includes open research questions, such as how to couple ad hoc file systems with the batch scheduling environment and how to schedule stage-in and stage-out processes of data between the storage backend and the ad hoc file systems. Also presented are strategies to build ad hoc file systems by using reusable components for networking and how to improve storage device compatibility. Various interfaces and semantics are presented, for example those used by the three ad hoc file systems BeeOND, GekkoFS, and BurstFS. Their presentation covers a range from file systems running in production to cutting-edge research focusing on reaching the performance limits of the underlying devices.Peer Reviewe